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ABSTRACT 

We study for the first time the environment of massive black hole (BH) seeds (~ 10 4-5 Mq) 
formed via the direct collapse of pristine gas clouds in massive haloes 10 7 Mq) at z > 6. 
Our model is based on the evolution of dark matter haloes within a cosmological TV -body 
simulation, combined with prescriptions for the formation of BH along with both Population 
III (Pop III) and Population II (Pop II) stars. We calculate the spatially-varying intensity of 
Lyman Werner (LW) radiation from stars and identify the massive pristine haloes in which it 
is high enough to shut down molecular hydrogen cooling. In contrast to previous BH seeding 
models with a spatially constant LW background, we find that the intensity of LW radiation 
due to local sources, Ji OC ai> can be up to ~ 10 6 times the spatially averaged background 
in the simulated volume and exceeds the critical value, J cr it, for the complete suppression 
of molecular cooling, in some cases by 4 orders of magnitude. Even after accounting for 
possible metal pollution in a halo from previous episodes of star formation, we find a steady 
rise in the formation rate of direct collapse (DC) BHs with decreasing redshift from 10~ 3 
Mpc~ 3 z _1 at z = 12 to 10~ 2 Mpc _3 z _1 at z = 6. The onset of Pop II star formation at 
z « 16 simultaneously marks the onset of the epoch of DCBH formation, as the increased 
level of LW radiation from Pop II stars is able to elevate the local levels of the LW intensity to 
•/local > Jait while Pop III stars fail to do so at any time. The number density of DCBHs is 
sensitive to the number of LW photons and can vary by over an order of magnitude at z = 7 
after accounting for reionisation feedback. Haloes hosting DCBHs are more clustered than 
similar massive counterparts that do not host DCBHs, especially at redshifts z >, 10. Also, 
the DCBHs that form at z > 10 are found to reside in highly clustered regions whereas the 
DCBHs formed around z ~ 6 are more common. We also show that planned surveys with 
James Webb Space Telescope should be able to detect the supermassive stellar precursors of 
DCBHs. 
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1 INTRODUCTION 

It is now an est ablished fact that galaxies host black holes (BH) 
at their centre s I Gebhardt et al. I l2000t iFerrarese & MerrittltZOoA 
iGultekin et alj|2009l) with BH masses ranging from 10 6 " 9 ' 5 M . 
The most massive BHs or supermassive black holes (SMBH) are 
believed to fue l quas a rs observed as early as z > 6 (see e.g . 



iFan et al.l [20031 . |2006| ; Iwillott et al.l [20031 : iMortlock et all 1201 lb 
This implies that the seeds of these SMBHs must have formed and 
grown to supermassive scales in the short time before the Universe 
was even one billion years old. It has also been suggested recently 
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dTreister et al.ll20lH Til here after) that there might be a popula- 
tion of obscured inter mediate mass black holes (IMBH s) at z > 6 
(however also see e.g. IWillottll201 J iFiore et aT who chal- 

lenge the claim). However, the origin of these SMBHs or IMBHs 
in the early Universe is still an open question. 

The most obvious way to make the SMBH seeds is from the 
stellar BHs in the early Universe. Detailed studies have shown 
that the first generation of stars (Pop III) form from metal-free 
gas, comprising mainl y of atomic and molecular hydrogen at earl 



times (see reviews by iBromm &~ Larson 2004]; ICiardi et 



n at early 
all 1200 ll 



and references therein). Pop III stars with masses in the range 
40 Mq < M» < 140 Mq and M» > 260 M^ colla pse into 
a black hole with M. =0.5 — 1 M* faeger et alj|2003l) and ac- 
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cretion of gas onto these stellar BHs offers a natural way to grow 
SMBHs, given their abundance and early formation times. 

This scenario however, has been challenged given that Pop III 
remnant BHs may not constantly accrete at or near the Edding- 
ton limit, which is likely required for 100 Mq seed black holes 
to reach a mass of 10 9 Mq b y z ~ 6. Both the radiation from 
the Pop III progenitor sta r (e.g. lYoshidalbood Ijohnson & Bromrnl 
l2007l : lAlvarez et al l2009h and the radiation emitted in the accretion 
process itself (e.g. lMilosavljevic etaDl2009l : IPark & Rico7till201ll : 
Il1|201 IJ), result in feedback which might slow down gas accretion. 
The constant availability of gas in the halo during the accretion 
period would also require the haloes to grow rapidly via mergers 
since episodes of star formation and feedback fr om supernovae ca n 
deplete the gas in such primordial haloes (e.g. iMori et al.ll2002h . 
On the other hand, a scenario where the accretion must be super- 
Eddington for a short peri od of time has been prop osed in order to 
allow fast BH growth (e.g. lVolonteri & Reesl200j) . which could be 
a result of the inefficient radiative losses due to the trapping of pho- 
tons i n the accretion disc (see e . g . IB egelmanll 1 9781 ; IWvithe & LoeTJ 
l20~il . 

Another possibility of growing stellar black holes is via merg- 
ers of haloes hosting either stars or BHs. A dense cluster or group 
of stars provides conditions under w hich frequent merg ers can oc- 
cur, leading to a runaway collapse dZwart et al .11 19991) that result 
in BHs with masses of around 10 J Mq. Mergers of Pop III seed 
BHs at hi gh redshifts can also, in principle, build up supermas- 
sive BHs fanaka & Haimanll 20091) . although slingshot effects and 
merger ti me scales pose problems for this scenario (see e.g. the re- 
views by|Natarajan 2 Tit I Volonterill201fjl and references therein). 

An alternative scenario is to make seed BHs with an ini- 
tial M. — 10 4 ~ 5 Mq via the direct collapse of pristine gas 
in haloes with T vir > 10 4 K (see e.g . lEisenstein & Loebll 1993; 



Oh & Haiman 2002 ; Bromm & Loe 



120031 : 



Koushiappas et al. 



2004; Lodato & Natarajan 2006). The key idea is to keep the haloes 



free of molecular hydrogen so that the gas collapses isothermally 
only via atomic hydrogen. For the gas collapse to proceed without 
fragmenting into stars, it also has to redistribute its angular momen- 
tum and various processes have been suggested in order to allow 
this, as explained below. 

In low spin haloes the gas settles down into a disc where the 
angular momentum can then be redistributed via gravitational insta- 
bilities, hence keeping the Toomre parameter c lose to unity and pre- 
ventin g the disc from further fragmentation dLodato & Nataraiar] 
120061 LN06 hereafter). The central core of M = 10 4 " 5 Mq, fed 
by the streams resulting from the non-axi symmetric disc instabil- 
ities, ultimately collapses into a BH with a similar mass. An im- 
portant feature of LN06 is that they explicitly link the dark matter 
halo properties, like spin and virial temperature, to the properties 
of the BH seed. Their model predicts the required ratio of the gas 
temperature to the virial temperature and the maximum halo spin 
which determines the final mass of the BH seed. 

The redistribution of angular mome ntum can occur via the 
'bars-within-bars' scenario as explored by iBegelman etail J2006h 
where the gas collapses into a dense self-gravitating core sur- 
rounded by an envelope supported by radiation pressure. The 
gas finally cools and collapses catastrophically via neutrino emis- 
sion into a central BH with an intermediate stage of a quasi-star 
jBegelman e"tai]|2008l). 



Spaans & Silk! d2006t ) showed that if the collapse of gas (com- 



of the medium. The time required for the Lyman-alpha photons to 
escape the medium becomes larger than the free fall time of the 
gas which prevents the gas from cooling and forming Pop III stars. 
Hence, the collapse can result in a massive BH which is of the order 
of 3 — 20 % of the total baryonic mas s of such haloes. 

Also, iRegan & Haetmeltl ( 120091) explored the gas collapse in 
rare atomic cooling haloes which could in principle host a DCBH 
in cosmological hydrodynamic simulation. They find cases where 
the inflow rates are high enough (> 1 Mgyr -1 ) to allow for the 
formation of massive BH seeds. 

All these scenarios end in a direct collapse black hole (DCBH) 
with M. ~ 10 4 ~ 6 Mq. Another alternative scenario includes the 
formation of a supermassive star (SM S) in an intermed iate step on 
the way to the formation of a DCBH l lBegelman|[2oTo]) . For this to 
occur the gas does not only need to be free of H2 and metals but the 
accretion rate onto the SM S needs to be high enough to allow th e 



•> needs to be nigh enough to allow the 
teegelmanll2010l ; Ijohnson et ai]|2012l) . 



prised of atomic H) in such haloes proceeds via an equation of state 
with a polytropic index larger than unity, Lyman-alpha photons can 
get trapped in highly dense regions owing to the large optical depth 



rapid growth to 10 4 ~ 6 M Q 

Although these scenarios take place in haloes with T v j r > 
10 4 K, which are mostly composed of atomic hydrogen, molec- 
ular hydrogen can form in these haloes when the densities are 
high enough to allow three-body hydrogen interactions. Such high 
particle densities are found at the halo centre and during the 
end stages of gas collapse. Hence these scenarios require a criti- 
cal level of H2 photo-dissociating Lyman Werner (LW) radiation 
(hv — 11.2 — 13.6 eV) in order to keep the abundance of H2 
molecules very low, as otherwise H2 cooling will lower tempera- 
tures to T ~ 200 K, thereby reducing the Jeans mass and leading 
to fragmentation of the gas cloud, which would ultimately result in 
star formation instead of a central BH seed. 

The main challenge in all the above DCBH formation scenar- 
ios is to reach the critical level of Lyman Werner radiation required 
to dissociate H2 molecules in the halo. Typical levels of a smooth 
uniform LW background, Jb g , range from 0.001 — 0.1 (where J 
is expressed in units of 10 -21 erg s~ 1 cm~ 2 Hz ~ 1 sr~ 1 ) and de- 
pends on the stellar density at a given redshift dGreif & Bromml 
1 20061) . whereas the critical value, J cr it, required for di rect collapse 
is ~ 30 (from Pop II) and ~ 1000 (from Pop III) JShang et alj 
l2010l ; IWolcott-Green et al.l201ll CS10 and WG1 1 herafter). It has 
been argued that a halo can be exposed to a radiatio n level higher 
than J cr it if it lives close to a star forming region dDiikstra et alj 
I2008L D08 hereafter). They use an analytical approach employing 
Poisson statistics and extended Press-Schetcher mass functions to 
model their halo distribution which accounts for clustering of the 
DM haloes and the spatial distribution of LW sources. 

Previous studies of DCBH forma tion have either assumed 
a spatially con stant LW background dRegan & Haetmeltl [20091 : 
IPetri et aT]|2012l) or a spatially varying LW background using ana- 
lytical prescriptions for clustering of sources D08. The latter study 
showed that the clustering of sources plays a crucial role in elevat- 
ing the levels of LW radiation above the critical value required for 
DCBH formation. While it is important to model the clustering of 
sources properly, it is also crucial to know whether a halo, which is 
exposed to the critical level of LW radiation, had previous episodes 
of star formation which enriched th e gas in the halo with metals. 
In contrast to D08, |Petri et al.1 d2012l) attempted to model the merg- 
ing histories of haloes using Monte-Carlo merger trees however, 
they did not account for the self consistent build up of the spatially 
varying LW radiation field. 

Due to the importance of LW feedback at high redshifts, 
some recent studies have exp lored the effects of LW radiation on 
early str ucture formation (e.g. Kuhlen et al ]|201ll). Pop III star for- 
mation dSafranek-Shrader et al.ll2012l IWhalen et alj|2008t iRicottj 
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20081). the evo lution of pair instability supemovae JWisell2012l ; 



Hummel et a l. 2011) and also on th e formation of SMBH seeds by 
Pop III stars dDevecchi et alj2012T) . 

In this paper we simultaneously follow the build up of the spa- 
tially varying LW radiation field as well as track the enrichment 
histories of dark matter (DM) haloes in a cosmological DM only, 
TV-body simulation using a semi-analytical model (SAM). We in- 
vestigate the conditions under which the LW intensity seen by an 
individual halo will reach a value £ J cr jt and we describe the re- 
sulting consequences for the formation of seed BHs via direct col- 
lapse. The aim of our work is to determine the plausibility of the 
existence of DCBH sites and probe the clustering features of such 
haloes. 

This paper is organised as follows. We describe the simulation 
and our model in the next section (Sec. [2} followed by which the 
results of our work are presented in Sec. [3] The observability of 
the supermassive stellar seeds of DCBHs by the James Webb Space 
Telescope (JWST) is discussed in Sec. [4] Finally we present the 
summary and discussion of our work in Sec. [5] 



2 METHODOLOGY 

In the following section(s) we describe our SAM that models the 
build up of the LW radiation field on top of our TV-body simulation. 
We model both Pop III and Pop II star formation and include a pre- 
scription for the evolution of the star forming and non-star-forming 
gas within an individual halo. This allows us to track the star for- 
mation histories of the haloes and account for the LW photon travel 
times which is needed in order to self consistently model the global 
and spatial level of the LW radiation at each point in our box. 



2.1 The TV-body simulation 

We use a DM only TV-body simulation with 768 3 particles in 
a 3.4 Mpcfo~ co-moving periodic box using the GADGET code 
dSpringel et al.l200ll : ISpringell2005l) . We assume a ACDM cosmol- 
ogy with n = 0.265, Q b = 0.044, fi A = 0.735, h = 0.71 and 
as = 0.801 consistent with the WMAP 7 results dKomatsu et alj 
l201lh . The resulting individual DM particle mass is 6500 M^/i" 1 
Merg er trees are constructed on the SUBFIND output dSpringel et alj 
120011) using the same method as in lSpringel et alj ( 120051) . Informa- 
tion on each subhalo includes its mass as assigned by SUBFIND, 
along with its host friends-of-friends (FoF) mass. The smallest re- 
solved DM halo contains at least 20 particles, which corresponds to 
1.3 x 10 5 Mq/i _1 . We run the simulation down to 2 = 6 and the 



snapsh ots are taken a few tens of Myr apart. As in ISpringel et alj 
J2005h and lCroton et al. U2006I) , our merger trees are based on sub- 
haloes. Note that we shall use the term halo instead of subhalo in 
the remainder of this work for the sake of simplicity. 

At a given snapshot, we label the haloes as minihaloes and 
massive haloes if their virial temperature is 20007^ ^ T v i r < 
10 4 K and T v i r ^ 10 4 K, respectively. Also, J or the combina- 
tion of the variable with any superscript/subscript explicitly implies 
Jlw in units of 10 -21 erg _1 s _1 cm _2 Hz _1 sr _1 unless specified 
otherwise. 

We define the infall mass, Minfail, of the halo as its mass at 
the last snapshot where it was the most massive subhalo within its 
FoF halo. We did this by tracking the halo's main progenitor branch 
back in time. The infall redshift is defined as the redshift when the 
infall mass was found. 

We use the relations from iBarkana &Loebl fcOOll) for the 



virial temperature, virial radius, T? v ir, and circular velocity, V c , of 
a halo 



T vir = 1.98 x 10' 



4 / M 
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(3) 



where fi is the mean molecular weight (1.22 for neutral primordial 
gas), fio is the matter density of the Universe at z = 0, Sl m {z) 
is the matter density of the Universe as a function of redshift and 
A c is the collapse over-density and z denotes the infall redshift as 
computed from our trees. 



2.2 Star formation 

In order for the first star to form out of the gas in a virialised pris- 
tine halo, the cooling time, t coo i, for the gas must be less than the 
Hubble time, tHubbic- The primordial gas mostly comprises of ei- 
ther atomic or molecular hydrogen and the cooling time depends 
on their respective cooling functions. Atomic hydrogen cooling is 
effective at T > 10 4 K whereas molecular cooling can operate at 
lower temperatures. In our model, since we probe the universe a t 
z ^ 30, we use the results from the study bv lTegmark et al.ldl997l) 
which showed that the critical fraction of H2 molecules required 
in order to satisfy the condition t coo i < tHubble is found in haloes 
with Tvi r ~ 2000 K at 2 = 25. Hence, the first star to form from 
a pristine gas cloud would be a Pop III star forming in a minihalo. 
The metals ejected from the first Pop III star would be enough to 
pollute th e gas and Pop II st ars could form subsequently in the same 
halo (e.g. lMaio et"ai1l2010l) . We discuss the Pop III and Pop II star 
formation in more detail in the following sections. 

As explained above, since it is critical to resolve minihaloes 
of mass ~ 10 5 ~ 7 Mq, this requirement limits the volume that we 
can probe with sufficient resolution in our study. We plot the mass 
functions of the FoF and subhaloes in our work at 2 = 6 in the 
Appendix. 



2.2.1 Pop III stars 

In our model, we allow a single episode of instantaneous Pop 
III star formation in pristine haloes with with T v i r ^ 2000 K 
dTegmark et al.ll997l : lM~aio et al.l20ldl) . Here, we consider a halo to 
be pristine if none of its progenitors have hosted a star in the past. 
In addition, our implementation of LW feedback, as explained in 
Sec. 12.31 regulates which pristine haloes form Pop III stars. Hence 
the non-Pop Ill-forming minihaloes can later to grow into pristine 

massive haloes. 

We assume a Salpeter IMF dSalpeterlll955l) with a mass range 
between 100 and 500 Mq and assume that on e Pop III star forms 
per minihalo (see e.g. iBromm & Larsonll2004l) . However, in mas- 
sive pristine haloes (T v i r Si 10 4 K) we form 10 stars following a 
Salpeter IMF, with mass cut offs at 10 and 100 Mq (e.g. lGreif et alj 
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l2008t|johnson et al]2008tlWise & Abell2008l) . Our choice of IMFs 
and mass cut-offs in both minihaloes and massive haloes is primar- 
ily to maximise the LW output from the stars. Forming multiple 
lower mass stars as opposed to a single very massive star gives an 
upper limit to the amount of LW radiation that can be emitted from 
a massive pristine halo as, for instance, the number of LW photons 
produced by five 100 Mq stars is larger than for one 500 Mq star 
(see Sec. |2.4t . 

Since the form ation time for a Pop III star is few Myr (e.g. 
iBromm et al]2009h and our snapshots are ps 10 Myr apart, Pop III 
stars are assigned a time of birth and distributed uniformly in the 
time interval between two subsequent snapshots (see Appendix). 
The masses of individual stars within a pristine halo are generated 
randomly following the respective IMFs assumed. 



2.2.2 Pop II stars 

The second generation of stars, Pop II, is also exp ected to exist 
at high redshifts within metal-enriched regions (e.g. IWise & Abeil 
l200alGreifetai]|201Ch . These stars are metal rich as compared to 
Pop III but have metallicities much smaller than the solar metallac- 
ity, Zq. The metals ejected from Pop III stars p ollute the host and 
neighbouring ha loes via stellar and SN winds dMori et alJIiool 
iMaio et al]|201 ll) . Any further collapse of the polluted gas in the 
haloes would result in cooling to lower temperatures, thereby re- 
ducing the Jeans mass and forming metal-enriched sta r s with lower 
masse s than the Pop III stars (e.g. IClark et al.l [20081 : ISmith et alj 
l2009h . The critical metallicity at which the transition occurs from 
Pop I II to Pop II ranges from 10~ 4 to 10~ 6 Zq (e.g. lFrebel et alj 
|2007|) . For simplicity, we consider a halo that has hosted a Pop III 
star (or merged with a halo hosting or having hosted a Pop III or 
Pop II star) pollute d with metals and a possible site for Pop II star 
formation (see e.g. |johnsonll2010l) . 

Since metals are the coolants required for making Pop II stars, 
we assume that a large enough potential well would be required to 
constrain the metals ejected from Pop III SNe and additionally add 
a constrain by setting t he threshold halo masaj for Pop II star for - 
mation to 10 8 M (e.g. lKitavama et alj2004l : lwhalen et alj2008h . 
Within these candidate haloes, we assume that the baryons can exist 
in either of the three phases i.e. non-star-forming gas, star-forming 
gas or stars. Below we describe the transition between these phases 
which ultimately regulates the Pop II star formation in a halo. 

• Non-star-forming gas phase: We assume in our model that 
once a DM halo crosses our resolution limit of 20 particles, it is ini- 
tially comprised of non-star-forming gas, A/hot = /bAfoM, where 
/b = 0.16 is the universal baryon fraction and Mum is the halo's 
current DM mass. 

While the DM halo grows between two snapshots, we add non- 
star-forming gas to the halo by calculating the accretion rate, M aC c, 
defined as 



AL 



JbAMoM — M* jP — Mout.p 
At 



(4) 



where AA/dm is amount by which the DM halo grows between 
two snapshots which are At apart, M», p and M out , p represent 
the total stellar mass and net mass lost in previous SN outflows 
summed over the incoming merging haloes respectively. 



1 See Section[5]for a critical view on how the choice of the Pop II threshold 
mass sets the clock for DCBH formation. 



• Star forming gas phase: In order for the gas to form stars, it 
must cool and collapse within the halo. We model the transition 
from the non-star-forming gas phase to the star-forming gas phase, 
-Mcdd, by allowing A/hot to collpase over the dynamical time of 
the halo, fd yn = ■% it - This estimate is justified by the fact that at 
such high redshifts, the radiative-cooling time is shorter than the 
dynamical time of the halo. 

• Star formation law: We then model the P op II star formation 
via a Kennicutt-type relation dKennicutl 19981) 

a 



AL 



ii 



-M, 



cold 



(5) 



O.ltdyn 

where a is the star formation efficiency (SFE). The factor 
O.ltdyn, which is the s t ar formation tim e scale, is motivated by 
iKauffmann et al] Jl999l) ; lMo etai] 1 19981) . 

Local observations indicate an a ~ 0.2, however, at this 
stage it is not clear if this also holds at high redshifts (z > 6), 
dKhochfar&Silkll201lh . We therefore treat a as a free parameter 
and normalise our model to the observations of the cosmic SFRD 
at z £ 6. 

• Outflows: In addition to star formation, we also consider the 
SN feedback processes in a star forming halo. We model the out- 
flow rate of gas from a Pop II star forming halo via the relation 



M out = 7 M,, n 
where 



7 = 



Vc 
Vout 



(6) 



(7) 



The functional form of 7 is taken from ICole et al l d2000h . We 
normalise the parameters in Eq. [7] to the results of the high reso- 
lution hydrodynamical simulations of the high redshift Universe 
(Dalla Vecchia and Khochfar 2012, in prep) and for the halo mass 
range considered in this work. This yields an outflow velocity 
Vout = 110 kms -1 and /3 = —1.74 resulting in typical values 
of 7 ~ 20. We assume that the outflows are generated in the 
star-forming gas phase and hence M out is subtracted from A/ co id ■ 

• Implementation: Each time interval between two consecutive 
snapshots, At, is divided into 100 smaller intervals and the follow- 
ing set of coupled differential equations (along with Eqs.[4][5]and 
[SJ for the individual baryonic components are numerically solved 
over the small time steps: 

M cold = ^-AL, u -A4 out , 



Mw = 



tdyn 

Mhot 

tdyn 



(8) 



(9) 



2.3 Impact of LW radiation on star formation and direct 
collapse 

Once the first generation of stars form in the Universe, the effects of 
LW radiation become important fo r subsequent star formation (e.g. 
Haiman et ^l200d : [Omukaill200lh . Even a small, uniform Jlw ~ 
0.01 from these stars can affect Pop III star forming minihaloes by 
dissociating a fraction of the H2 molecules and preventing the gas 
from cooling and collapsing jMachacek et alJl200l1: lYoshida et alj 
120031 : IWise & Abe3l2007l : b'Shea & Normanll2008l) . The amount 
of H2 molecules that can be dissociated depends directly on the 
LW background it is exposed to. 



© 0000 RAS, MNRAS OOCUJfJJ] 



Ubiquitous seeding ofSMBHs by direct collapse 5 



10 7 



10 e 



0.001 



□ □ 



Pop III forming Minihaloes 



0.010 



Non star forming Minihaloes 



0.100 



1.000 



10 4 ' 



10 4 



2L 10 42 



10 4, 



Cal.from Starburst99 
M. = 10 e M s 
Z = 0.001 Z s 
Salpeter [1,100] M e 



200 



400 600 
Age [Myr] 



800 



1000 



Figure 1. The fit used in our model based on the lO'Shea & NormanlfcOOlil) , 
Fig. 3(c) (square-symbols) which determines the minimum mass of a pris- 
tine minihalo that can host a Pop III star for a given level of external LW 
radiation it is exposed to. For a metal-free minihalo exposed to a given Jlw. 
if its mass lies above the line it is considered Pop III star forming. The solid 
line represents Eq. llQI 

The minimum mass, M cr i t , of a pristine halo in which the gas 
is able to cool, collapse and form Pop III stars in the presence of a 
given external LW radiation intensity can be approximated by 



M^t = ip 1.25 x 10 5 + 8.7 x 10' 



/ Jlw \ ' 



(10) 



where the expression within the brac kets is the functional fit to the 
numerical simulations carried out by iMachacek etail d200ll) . The 
correction fact or ip has been set to 4 follo wing the higher resolution 
simulations of lO' Shea & Normanl ( 120081) . as shown in Figure [TJ 

One might argue from Fig.[TJthat Jlw = 1 is sufficient to set 
the threshold mass to 10 7 Mq, which is the mass beyond which 
pristine haloes can cool via atomic hydrogen and hence, direct col- 
lapse should ensue. However, detailed simulations by CS10 and 
WG1 1 show that only the H2 molecules in the outer regions of such 
halos are dissociated and a considerable fraction of molecular hy- 
drogen (~ 10 -3 ) still exists in the central region of the halos in the 
presence of such low levels of Jlw - In order to prevent star forma- 
tion in the central parsec region of the halo, it is essential to bring 
down the H2 fraction in the gas to 10~ 8 which can be achieved by a 
J™ = 1000 from Pop III stars (WG1 1) or J° it = 30 - 300 from 
Pop II stars (CS10). The difference in the values of J cr i t for Pop 
III and Pop II stars is due to the difference in the spectral shapes 
of the two stellar populations. As shown by CS10, the lower value 
of J^ it can be partly attributed to the fact that the H~ dissociation 
rate from Pop II stars is ~ 10, 000 times larger than that from Pop 
III stars, due to the softer shape of the Pop II spectrum at 0.76 eV. 
Since, H~ is a precursor to H2 formation, destruction of is crit- 
ical as it results in a lower rate of H2 production^ 

Thus, if a metal free halo with T v i r > 10 4 K is exposed to a 



2 H is is dissociated by the following photoreaction: 
H~ + hv -> H + e~ 



(11) 



The dissociation rate can be written as &28 = 10 — 10 s~ 1 a Jlw- Here, 
Q III = 0.1 for Pop III stars and ot\i = 2000 for Pop II stars. Since H - 
can lead to H2 formation, this reaction is of prime importance in order to 
keep the gas at a low H2 fraction. 



Figure 2. Emission in the LW band from a Pop II population as a func- 
tion of its age comp uted using the data from STARBURST99 catalogue 



non ot its a ge computed usin 
iLeitherer et al.l 19991 Fig. 7e). 



critical level of LW radiation, a direct collapse can ensue. In this 
scenario, the cooling is suppressed and the gas stays at ~ 8000A'. 
Due to the large Jeans mass and high accretion rates that these high 
temperatur es imply, a SMS forms and subsequently co llapses into 
a BH (e.g. iBromm & Loeblliool : Ijohnson et all 120121) . The cen- 
tral BH then continues to accrete, embedded in an envelope of gas, 
at super-Eddington efficiencies reaching a quasi-star state and col - 
lapsing into a M. = 10 4 " 5 M black hole faegelman et al.l200ot) . 

To summarize, only metal-free minihaloes with masses larger 
than Merit (see Eq. UOt are considered to be Pop III star-forming. 
Also, if these minihaloes are exposed to a Jlw 1 then they 
are considered to be non Pop III star-forming. The metal-free mas- 
sive haloes still make Pop III stars irrespective of Jlw, given that 
Jlw < Jcrit otherwise they can be considered to be DC candidates. 
Also, the Pop II haloes are unaffected by any value of Jlw since 
they are polluted by previous generations of stars, have Mi n faii > 
10 8 Mq and the coolants are metals. Therefore, the Pop II star- 
forming criteria for these haloes is that Mi n f a ii > 10 s M Q and that 
the halo has hosted stars previously or has undergone a merger with 
a previously star-forming halo. 

Hence, for pristine minihaloes, i.e. haloes with masses in the 
range corresponding to 2000 ^ T v i r < 10 4 K, 



M > M cr it 

Jlw < 1 



Pop III 



The other pristine minihaloes that do not satisfy the above condi- 
tions can not form Pop III stars. 

For pristine massive haloes i.e. haloes with masses corre- 
sponding tO Tvir > 10 4 K, 

Jlw < Jcrit } Pop III 

Jlw ;S Jcrit } DCBH 



2.4 Jlw calculation 

We describe our calculation of the mean and local LW intensities 
from both the Pop III and Pop II stellar populations in this section. 
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Table 1. Functional fits used to co mpute the radiat ion output and age of 
Pop III stars. Values are taken from[Schaerer (2Qo2) for a zero metallicity, 
no mass loss case. 



X 


a 


ai 


a 2 


a 3 


age (Myr) 


9.785 


-3.759 


1.413 


-0.186 


Qh 2 


44.03 


4.59 


-0.77 




Qh 


43.61 


4.90 


-0.89 





2.4.1 Mean Jlw calculation 

The first stellar populations in the Universe mark the onset of the 
ultra-violet (UV) background which has a negative effect on star 
formation as described in the previous secti ons. The LW photo n 
horizon is larger than our box size (~10 Mpc. lHaiman et alj2000h . 
Therefore the contribution to the background must come also from 
outside our simulated volume. In order to account for this, we as- 
sume that the SFRD in our volume is representative of a larger cos- 
mological volume. The mean LW background in our volume, is 
then assumed to exist everywhere in the Universe and is assumed 
to be the minimum level of LW radiation that a halo is exposed to 
at any given redshift. It can be computed following the formulae in 
iGreif & Bromrril d2006l) : 

J™ - /osct--— »7lwP» x (l + Z Y 



J, 



f, 



47rraH 
he 



(13) 



' 47T771H 

where / csc is the escape fraction of LW photons from the halo, 
pi 11 , pj, 1 denote the comoving density of Pop III and Pop II stars 
respectively at the given redshift z, ??lw is the number of LW pho- 
tons per stellar baryon (t/lw — 10 4 a nd t?lw = 4 x 10 3 fo r the 
assumed IMFs in our study and as in IGreif & Brommll2006l) and 
h, c, rriH are the Planck's constant, speed of light and mass of a 
hydrogen atom respectively. In our model, the parameters p™, p 1 } 
are computed by checking if a star or stellar population is active 
at the current snapshot. Each stellar source in our model is given a 
time of birth, which is the epoch at which the star is formed and a 
lifetime depending upon the mass (see Appendix for more details). 



2.4.2 Spatial variation of Jlw 

It is important to note that Eqs. 1121 and 1131 are valid for a mean, 
uniform LW background. However, it is possible that a halo would 
have some stellar sources in neighbouring haloes which would pro- 
duce levels of Jlw higher than the global mean value, which would 
depend on the clustering scale of haloes dAhn et al]|2009i KA09 
hereafter). 

In order to calculate the effects of a spatially varying Lyman 
Werner specific intensity from individual Pop III stars, we write 



jIII /csc fo^avg 



Qlw 



TY Az/LW 47TtP 



(14) 



were /t^avg is the average energy of a photon emitted from a Pop 
III star in the LW band, A^lw is the difference in the maximum 
and minimum value of the LW frequency range, d is the luminosity 
distance, and Qlw (expressed as; Qlw = Qh 2 — Qh) is the num- 
ber of photons produced per second in the LW energy range. The 
factor of ty in Eg. 1 141 T ^77^ ) arises from the conversion of the flux 
into specific i ntensity, assuming that each Pop III star is a uniform 
bright sphere dRvbicki & Lightmanlll986l ). The specific values that 
we use for these parameters in the case of Pop III stars have been 



2 
1 


-1 
-2 
_ -3 
-4 
-5 

-3 -2 -1 

log SFR [M e yr 1 ] 

Figure 3. Comparison between the SFR-function as computed from S12 
and our work. The "I'sfr, bounds at z = 6 and 7 as computed from S12 
are represented by the red and blue regions whereas the data from our sim- 
ulation is marked as the red (z = 6) and blue (z = 7) filled circles. The 
current observational surveys are able to probe t he SFRs rightwa rds of the 
vertical lines denoting SFR ~ IMq yr -1 fe.g. lSmitet all2012h . 




computed using the functional fits from lSchaererl J2OO2I) where they 
track the evolution of stars with different masses and metallicities 
in their models. They express a given parameter X (see Table[T} as 
a function of the stellar mass M*, as follows: 

log(X) = ao + aim + 02m 2 + 03m 3 , (15) 

where m = log(-£M. 

For the contribution of LW photons from Pop II stellar 
sources, since we only form a total mass in Pop II stars (Af»,n, 
see section |2. 2. 2t . we calculate Pop II properties using the data^f] 
obtained from the STARBURST99 catalogue. We integrated the 
curve(s), in the LW range to obtain a function E which is the en- 
ergy per unit time (in units of erg sec" 1 ) emitted by a 10 6 Mq Pop 
II stellar population as a function of the age, as shown in Fig. [2] We 
then calculate 



T n - 

"■'local 



E M fi 



TY A^LW 47TC? 2 



(16) 



where Me.*,n is the mass of the stellar population normalised to 
10 6 M Q . Similar to Eq. [H the factor of ty in Eq. [T6l 
arises from the conversion of the flux into specific intensity, as- 
suming that the Pop II stella r population is a uniform bright sphere 
dRvbicki & Light man 1986). We form Pop II stars following the 
prescriptions described in Sec. 12.2.21 As discussed in Section l2.4. II 
we assume each halo is exposed to a minimum level of J = Jb g 
calculated in each of our runs at each snapshot. We then add up 
the LW contribution from each star on top of the background to get 
the total value of the LW radiation that a halo is exposed to. This 
slightly overestimates the LW contribution by a factor of less than 
a few percent. However as we will show later this will not impact 
our results. To summarise, we have 



J 



•A 



J, 



(17) 



3 The data from STARBURST99, Fig. 7e assumes a Salpeter IMF with a 
mass cut off at 1, 100 Mq, instantaneous star formation, total stellar mass 
= 10 6 Mq, metallicity of Z = 0.001 and no nebular emission. These 
parameters are the closest to a Pop II stellar population. 
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Observations 
Poplll(esc0.5) 
Pop II (escO.5) 

Pop III (other cases) 
Pop II (other cases) 




Table 2. Summary of cases considered in our work. The fiducial case is 
marked in bold. 
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Figure 4. SFRD computed using the methods described in Sec. |2]for all 
the cases in our work (see Table. |2j- Solid, dotted, dashed, dash-dot-dash 
and dash-double-dot-dash represent the cases esc0.5, escl.0, esc().5HSFE, 
escO.l and esc0.5reion respectively. The green square, upright triangle, 
inverted triangle and cir c les represent observat ional data from iHopkinsI 
<2004 . lMannuccil l2007t) . iBouwens et alj liooij) and lLaporte et al.N2012h 
respectively. 



.II .II 

"^he- ~r t/local ; 



J — Jbg 

T - T UI _L 7 

Jtotal — J + J 



(18) 
(19) 



The quantity Jtotai is only used for determining if the pristine 
minihaloes can host Pop III stars (see Fig. QJ. In their work, 
lO'Shea & Norman] d2008h analysed the gas collapse within haloes 
in the presence of a Jlw flux. The photons could be coming from 
Pop II, Pop III or both as long as the photons are in the correct en- 
ergy band, hence Eq.[l9]is valid for analysing pristine minihaloes 
for Pop III star formation. 

On the other hand, due the importance of the spectrum at lower 
energies for the dissociation of H~, the quantities J m and J u are 
used to determine if the gas in the halo can undergo DC by com- 
paring the values to J™ t and Jlj it respectively. 



2.5 Escape fraction of LW radiation and reionization 
feedback 



Recent studies 
IPaardekooper et all 



i Wise & CerJ 120091 : 
shown 



201 ll) 



have 



lYaiima et alj 1201 ll ; 
that that the escape 



fraction for UV photons could vary with the parent halo mass, 
however, the precise values of LW escape fractions from haloes is 
still unclear. One might argue that once a pristine halo has hosted 
a Pop III star (or even a Pop II stellar cluster), most of the H2 
is depleted in the halo and the LW photons should, in principle, 
escape the halo unobstructed, implying / csc ,haio — 1-0 (KA09). 
However, the stars (which form in dense environments within 
the halo) are expected to be surrounded by molecul ar hydrogen, 
hence implying a /esc.haio < 1-0. Previous studies , jRicotti et alj 
l200ll : iKitavama et alj[2004ri have found that the minimum escape 
fraction for LW photons can be 0.1 but can also reach values of > 
0.8 in minihaloes. 

In addition to the / ea c,haio, the optical depth, tlw, of the inter- 
galactic medium (IGM) would also i mpact the number o f LW pho- 
tons reaching a neighbouring halo. ICiardi et alj (120000 find that, 
typically, tlw ~ 3, and including this in our calculations would 
imply an additional factor of e~ rLW in Eqs. [HI [HI [HI and [76l 



Case 


/esc 


Feedback 


a 


escO. 1 


0.1 


LW 


0.005 


escl.0 


1.0 


LW 


0.005 


escO.S 


0.5 


LW 


0.005 


esc0.5HSFE 


0.5 


LW 


0.1 


esc0.5Reion 


0.5 


LW + reionisation 


0.005 



Note that effectively, the / csc used in our work can be viewed 
as a degenerate combination of an escape fraction of LW pho- 
tons from the halo, / esc ,halo> arid the optical depth of the IGM i.e. 

fcsc = ,/csc,halo X € Lw . 

Given the uncertainty in / OS c,haio and tlw. we chose three 
cases to bracket the range of possibilities: / osc = 0.1,0.5,1.0. 
We also ran a case for our model in which we implemented an 
additional reionisation feedback from hydrogen-ionising photons 
by setting a circular velocity threshold of 20 kms -1 for all the 
haloes with T v i r > 1 4 K at 6 < z < 1 0. This choice is mo- 
tivated by the work of iDiikstra et alj (|2004j) where they study the 
gas collapse in haloes under a photo-ionising flux and find that 
a halo must be above a certain mass threshold (characterised by 
circular velocity in their work) to allow fo r at least half of the 
gas to undergo collapse. O ther studies (e.g. lOkamoto etaHl2008l : 
llVIesinger & Diikstral2008l) have also looked into the feedback pro- 
cess and have found similar mass thresholds. 



2.6 Model normalisation 

For our fiducial case, esc0.5, we set / CS c=0.5 for the LW radiation, 
a=0.005 and implement LW feedback in the model. The choice of 
a is made in order to match the observations of the cosmic SFRD at 
z £ 6. We normalise our free model parameter for the st ar forma- 
tion efficiency against recent observations o f the SFRD l|Hopkinsl 
|2004 lMannuccill2007l ; IBouwens et "alll2008l : lLaporte et alj|2012h . 



Due to the sensitivity limits of present surveys, the range in star 
formation rates probed in our simulations is not observed. Thus we 
chose the fiducial value of a in our model to lie within the error lim- 
its of the extrapolated faint-end slope of the SFR-funct ion, <E>sfr, 
of star forming galaxies at z = 6 and 7 as shown in ISmit et all 
j2012l) (S12 hereafter). In Fig. [3] the red and blue regions enclose 
the limits on 'I'sfr at z — 6 and 7 respectively, constructed using 
the fit parameters provided in S12. The blue and red filled circles 
denote the data from our work which is in fair agreement with the 
expected values of $sfr. 

We vary / osc to 0.1, 1.0, keeping a constant at 0.005 and 
name the cases escO. 1 and esc0.5 respectively. The model labelled 
esc0.5Reion is where we also account for reionisation feedback ef- 
fects on top of the fiducial case. We also implemented a high star 
formation efficiency for Pop II stars by setting «=0.1 and label the 
model as esc0.5HSFE. All our cases are summarised in Table [2] 



3 RESULTS 

The evolution of the SFR density (SFRD) with redshift for all cases 
is plotted in Fig. [4] The green symbols represent the cosmic SFRD 
as inferred from observations. The SFRD we compute is within the 
observational constraints at z ~ 6. Pop II stars first appear in our 
box at z ~ 16. The red and blue solid lines in the plot represent our 
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Figure 5. The self consistent build up of Jb g , defined as the addition of the 
LW background radiation from both the stellar populations, plotted against 
redshift. In all the cases except esc0.5reion, higher Jb g leads to a higher 
number of total DCBHs due to the efficient LW feedback and the resulting 
higher J t otal- 

fiducial case of esc0.5 and the light blue and orange lines represent 
the other cases in our work. The Pop II SFR is roughly the same 
in all the cases except for esc0.5HSFE and esc0.5Reion. However, 
due to the different escape fractions assumed for the cases, the Pop 
III star formation varies over all redshifts. This is due to the fact 
that in our model, although the Pop III star formation is critically 
affected by the self consistent build up of the LW radiation, Pop II 
star formation is not. In general, the Pop III SFRD is inversely pro- 
portional to the number of LW photons produced, which is directly 
proportional to the escape fraction. This is illustrated by the higher 
level of the Pop III SFRD in the escO. 1 case (light blue, dash-dot- 
dash line) than all the others. 

After z ~ 16, the increase in the LW radiation due to the Pop 
II stars (see the following sections) is able to further suppress the 
Pop III star formation. Also, maximal suppression of Pop III star 
formation is observed in the esc0.5Reion case where the additional 
mass constraint of 14 = 20 km s~ between 6 < z < 11 prohibits 
the pristine minihaloes and massive haloes from making Pop III 
stars. The additional circular velocity threshold in the esc0.5Reion 
case also causes a drop in the Pop II SFRD at 2 < 11 which is 
again due to the fact that the halo mass corresponding to V c ^ 
20 kms" 1 , the assumed mass threshold for structure formation in 
esc0.5Reion, is slightly higher than 10 s Mq which is the mass 
threshold for Pop II star formation in all the other cases. 

3.1 The LW intensity 

We start by expressing the total mean LW intensity at a snapshot 
as the sum of the contribution from the Pop III and Pop II stellar 
sources. 

Jbg = Jbg + Jbg ■ (20) 

The evolution of Jb g for all our cases is plotted in Fig. [5] Note 
that in each case, Jb g scales as the product of the escape fraction 
and SFRD at a given redshift. This is illustrated by the fact that 
although the level of the Pop II SFRD for esc0.5 and escO. 1 is simi- 
lar, as seen in Fig.[4] the level of Jb g is lower for escO.l than esc0.5 
in Fig. [5] The higher level of Jb g for esc0.5HSFE than esc 1.0 for 
11 < 2 < 16 can be explained by the SFRD in the respective 



cases. The SFRD is considerably higher between 11 < z < 16 for 
esc0.5HSFE, however once the SFRD approaches that of esc 1.0 
(not visible in Fig. [4] as it is hidden by the red line), the escl.O 
case produces more LW photons and hence a higher value of Jb g 
is seen for escl.O at 2 < 11 as compared to esc0.5HSFE. The 
build up of Jbg is in sync with the the SFRD in each of the cases 
and hence, consistent with our implementation. In Fig. [6] we plot 
the individual backgrounds from both stellar populations and the 
maximum level of the LW flux seen in a pristine halo at each red- 
shift, for our fiducial case. The solid blue and red lines represent 
J™ and J^g respectively which add up to the solid black line, Jb g . 
The Jbg computed at each timestep is assumed to be the minimum 
level of LW radiation to which a halo is exposed at that timestep. 
As expected, the maximum local level of the LW flux for a stellar 
population is always higher than the background (and in some rare 
cases equal to the background level). At all redshifts (except two 
cases at z ~ 12, 30 where the J™ t is seen by minihaloes), Pop 
III stars produce a Jm<J™ t . On the other hand, Pop II stars are 
able to produce a Ju>jH it in at least one of the pristine haloes at 
all redshifts, which is shown by the the red triangles being above 
the red dashed line. The epoch of DCBH formation in each of our 
cases is only observed after the Pop II star formation kicks in. 

We plot the distribution of the local Jlw as seen by pristine 
haloes in Fig. [7] before and after the Pop II star formation begins 
at z ~ 16. We define / pr i s as the fraction of pristine haloes with 
T v i r > 2000 K exposed to a given Jlw- The red and blue solid 
histograms represent J 11 (Eq|18t and J m (Eq|17> respectively. In 
Fig. [7] we see that less than one percent of pristine haloes (which 
roughly translates into a fraction of 5 x 10~ 4 of the total number 
of haloes) see a J 11 > Jlj it whereas the Pop III LW flux is always 
subcritical even before the Pop II star formation begins. The low 
fraction of pristine haloes that are exposed to j!J it can be attributed 
to the rarity of the event where a pristine halo is clustered (hence 
close enough) to the neighbouring haloes hosting Pop II stars. The 
trend of the distribution function is similar to D08, where they plot 
the PDF of all the haloes exposed to varying levels of Jlw- The 
Fig- Ul further supports our result from Fig. |6]i.e. Pop III stars are 
always subcritical to DCBH formation and that the LW radiation 
required for DCBH formation is always produced by Pop II stars. 

The value of Jb g that we compute in our work for esc0.5 
and esc0.5Reio n is w ithin 5% of the previous estimates of 
iGreif & Bromml ( 120061) at 2 = 10, where they self-consistently 
study the impact of two types of stellar populations and their feed- 
back on star formation at 2^5. The value of J™ that we find 
for all our cases, de noted by the so l id line s at z Si 16 in Fig. [5] is 
also consistent with Ijohnson et all d2008l) . i.e. it does not exceed 
their value of the maximum level of the LW background expected 
from Pop III stars (~ 0.13 at 2 ~ 16). We also find a good agree- 
ment with lTrenti & Stiavellil J2009l) as our SFRD and J outputs re- 
semble their estimates, but it is difficult to draw exact comparisons 
as they used an analytical Press-Schetcher modelling and the pa- 
rameter choices of the studies differ considerably. In cases esc0.5 
and esc0.5Reion, we find Jb g ~ 1 at 2 ~ 10 (also see Fig. [5j 
whic h is very close to the expected value during the reionization 
erfl l lOmukaill200ll ; iBromm & Loebll2003l) . 



4 In a recent attempt to model the DCBH formation. iPetri et alJfeOl j) find 
a very high global LW flux of Jlw — 1000 at the epoch of reionisation and 
they argue for even higher levels in spatially clustered regions. Our self- 
consistent methods to calculate the background and spatial variation of the 
LW flux show that such a high background is difficult to achieve. 
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Figure 6. The mean and maximum level of LW radiation plotted for each redshift in our fiducial case. The red triangles (J 11 ) and blue crosses ( J m ) indicate 
the maximum value of LW radiation to which a pristine halo is exposed at each redshift in the simulated volume. The red and blue dashed lines represent J,5J; t 
and J*H t respectively. It is interesting to see that the maximum value of J In (blue crosses) falls short of J™ t (blue dashed line). However, in the case of 
Pop II sources, the maximum value of J 11 (red triangles) is several orders of magnitude higher than the JH it (red dashed line). Hence we deduce that in our 
simulation, Pop II sources are most likely the ones to produce the J cr i t required for direct collapse. Note that the spatial LW radiation is computed only in the 
pristine haloes with T vir > 2000 K. 



In accordance with KA09, we find that the local LW intensity, 
which can be orders of magnitude higher than the mean LW in- 
tensity, is observed in highly clustered regions. This becomes even 
more evident in Sec. 13.4. fl where we present the cross-correlation 
functions of DCBH haloes. Note that although KA09 carried out 
full radiative transfer cosmological simulations to model the spatial 
variation of the LW intensity, they lacked the resolution required to 
study the impact of LW radiation feedback on structure formation 
(namely star formation in minihaloes) at such high redshifts. 



3.2 Sources responsible for J > J cr it 

With the change in the relative fraction of Pop III to Pop II star for- 
mation, the relative contribution to the LW background undergoes 
a change as well. As seen in Fig.[5]|6]and[7] the main contribution 
within our model comes from Pop II stars. At almost all times the 
contribution from Pop III stars is subcrititcal for DCBH formation. 
The critical level of radiation required for direct collapse is always 
produced by Pop II stars. This can be attributed to the way in which 
Pop III and Pop II stars form. We checked the effect of allowing Pop 
III stars to form as single stars, in binary systems or in a group of 
10 stars (which maximises the LW flux from a halo) with varying 
IMFs, and in all the cases the total number of stars was insufficient 
in producing the critical LW flux in a neighbouring halo. 

Also, the short lifetime of a Pop III star poses a problem as 
they reach the end of their lifetime in a few 10 6 yr as compared to a 
Pop II stellar population which can actively contribute towards crit- 
ical levels of LW radiation for up to a few 10 7 yr (see Fig. [2}. Hence 
even if it could produce J cr it, a Pop III star is less likely to be near 
a massive pristine halo (and hence contribute towards J cr it) than 
a Pop II stellar population. The result is in good agreement with 



llnavoshi & Omukail JioTlh , where they also argue that for a stellar 
source to produce the J or it, it must be a Pop II/I star cluster or very 
top heavy Pop III galaxies. As per our current understanding, Pop 
III stars form in (at most) groups of a few with masses ~ few tens 
of solar masses (e.g. iGreif et aTll201ll) . The occurrence of Pop III 
galaxies at such high redshift is expected to be extre mely rare as 
metal pollution in Pop III hosting haloes is quite fast l lMaio et al] 
l201ll) and it is highly unlikely that a cluster of these short lived 
Pop I II stars could end up in a galaxy Jjohnson et al .l2008l ; |johnsor] 
|2010|) . If a massive, pristine halo is to undergo direct collapse, we 
would expect it to have an external close by neighbour hosting a 
Pop II stellar population giving rise to Jcrit- 

The variation in the spatial LW intensity we find is also con- 
sistent with KA09. Using a full radiative transfer prescription in a 
cosmological box, they found that LW radiation varies on the clus- 
tering scale of sources at high redshifts, but they did not resolve the 
Pop III star forming minihaloes important for such studies. 



3.3 Abundance and growth of DCBHs 

The rate at which DCBHs are found to be forming in our simula- 
tion volume is shown in Fig. [8] We find a steady rise in the DCBH 
formation rate density with decreasing redshift (in units of Mpc -3 ) 
which can be expressed as 



dN 
dz 



bi (1 + z) 



(21) 



where 61, 62 are the fit parameters for the DCBH formation rate in 
each of our cases as shown in Table|3] As we are neglecting po ssi- 
ble metal pollution from neighbouring halos dMaio et alj|201 ll) the 
formation rates are strict upper limits in each of the cases. The fact 
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Figure 7. The distribution of Ji oca i for all the pristine haloes at a given 
redshift for our fiducial case. The red and blue solid histograms (dashed 
lines) indicate Jtt (Jff 1 ) and Jm (Jjjj 1 *) respectively. The absence of the 
red histogram in the top panel is because Pop II star formation begins at 
z = 16 in our box. 



that we find a few DCBH candidates in our 3.4 Mpch^ 1 simula- 
tion volume implies that the conditions for a DCBH are achievable 
in the early Universe and many such intermediate mass BH (if not 
SMBH) should exist at high redshifts. 

The DCBH formation rate increases as a function of the num- 
ber of LW photons that are emitted from a star forming halo. In 
escO.l we find 13 times fewer DCBHs than in esc0.5, which is due 
to the lower escape fraction assumed for LW photons in the former 
case. As the escape fraction increases from 0.5 to 1.0, the DCBH 
formation rate increases considerably. In all the above cases, this 
can be explained by the effect arising from the change in J to tai 
(Eq. 1 1 9I> which is two fold 

• a higher J to tai implies that more minihaloes are prevented 
from Pop III SF due to efficient LW feedback which makes them 
available for DCBH formation at later times since they are not 
metal-enriched. 

• a higher Jtotai directly affects the efficiency of DCBH forma- 
tion as it easier to exceed J cr it- 

The lower formation rate of DCBHs in esc0.5HSFE than esc 1 .0 can 
be attributed to the lower level of J Dg in the former case at z < 11. 
Also, since the majority of DCBHs form at z < 11, the fits are 
dominated by the DCBH formation rate at later times and hence a 
lower slope for esc0.5HSFE than esc 1.0 is seen. 
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Figure 8. Formation rate density of black holes in all the models plotted 
against redshift. The line represents a fit (Eq |21l to the formation rate den- 
sity in each case. 



The esc0.5reion case produces an interesting outcome where 
we find only 4 DCBHs which is roughly 10 percent of the DCBHs 
produced in our fiducial case. Note that before the reionisation 
feedback kicks in at z — 11, both esc0.5 and esc0.5Reion have the 
same DCBH formation rate. One would expect photo-ionisation ef- 
fects and the photo-evaporation of pristine minihaloes in the early 
Universe, which is accounted for by our reionisation feedback 
model, to greatly reduce the number of haloes into which primor- 
dial gas can collapse at later times. We discuss this case in more 
detail in Sec. 13. 61 

To explore the impact of BH growth via accretion after their 
formation, we allow the BHs to grow via Eddington accretion using 
the relation 



M,(t) = M.,o exp / dd 



t 



450 Myr 



(22) 



where M. is the final mass of the black hole, M m> o is the ini- 
tial mass of the bl ack hole set to 10 4 Mq for a DCBH (e.g. 
Ijohnson" et al.ll2012l CS10), t is is the accretion time, e is the ra- 
diative efficiency and / c dd is the Eddington fraction. We explore 
the (/odd, e) parameter space by choosing e — 0.07, 0.1, 0.2 to ac- 
count for a range in radiative efficiencies for Eddington (/odd = 1), 
sub-Eddingtion (/ c dd < 1) and super-Eddington (/ c dd > 1) accre- 
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Table 3. Fit parameters lines in Fig.[8]which follow Eg 12 H and total DCBH 
number density at z = 6 in each of our cases 



Case 


bi 


f>2 


Total DCBH 

(Mpc~ 3 ) 


escO. 1 


0.035 


-1.69 


0.13 


esc0.5 


0.054 


-1.79 


0.518 


escl.O 


49 


-4.4 


2.09 


esc0.5HSFE 


0.48 


-2.43 


1.58 


esc0.5Reion 


0.0061 


-1.0 


0.035 



tion jjohnson et alj|201 ll ; IShapirol2005l) . Since there is a lot of am- 
biguity regarding the early regimes of BH accretion, we varied our 
Eddington accretion parameters (/odd = 0.4, 1.0, 1.5) to account 
for a range of possibilities in the overall accretion mode of the BH. 

For our fiducial case (esc0.5), Fig. |9] shows the DCBH mass 
function constructed for / c dd = 0.4, 1, 1.5 and e = 0.1 (top panel) 
and the cumulative mass density of these DCBHs at z — 6 (bot- 
tom panel) plotted for different choices of / e dd and e. It is clear 
from the top panel that the DCBHs can almost reach SMBH scales 
with Eddington accretion and quite easily attain a mass larger than 
10 9 M sun if we assume super Eddington accretion. A wide range 
of evolved BH mass densities is seen in the bottom panel. Each 
filled black circle in the bottom panel represents the mass density 
of newly formed DCBH at that redshift, assuming an initial DCBH 
mass of 10 4 Mq. The solid purple triangle in the bottom panel is 
the observational claim made by Til for the mass density of the 
IMBH at z ~ 8. Although we do not match Til's claim for esc0.5, 
we do so for esc0.5Reion as explained in section [3~6l 




log(M)[M e ] 




3.4 DCBH Host Haloes 

In the following sections, we explore the regions where we find the 
conditions for direct collapse and the histories of the DCBH host 
haloes. 



3. 4. 1 Environment 

In order to understand the environmental differences between the 
haloes that host DCBHs and the ones that do not, we construct 
cross-correlation functions for the distribution of sources around 
them. In our case, the cross-correlation function is the excess prob- 
ability of encountering a source in a given distribution of sources 
around a halo as compared to a uniform distribution. Consequently, 
we define the cross correlation function as 



£(d) = 



DD(d) 
RR(d) 



(23) 



where DD(d) represents the data-data pair counts at a given dis- 
tance d, constructed from our model. The data-data pairs in our 
work refer to the halo-source pairs. We fix the halo and loop over all 
the qualifying sources thereby computing the physical distances. 
This is done for all the halo-sources pairs in a given redshift range. 

RR(d) represents the random-random pair counts constructed 
from a un iform dis t ributio n of sources around a random halo. 
Similar to iLi et alj d2012l) . we use the formulation of £(d) to 
qualitatively compare the small scale clustering properties of 



5 The qualifying sources for a given halo refer to the ones that satisfy the 
conditions described in section lA2l of the Appendix. 



Figure 9. DCBH mass function and cumulative mass density for the fiducial 
case (esc0.5). Red (/ cdd = 1.5), blue (/ cdd = 1) and green (/ odd = 0.4) 
represent the super-Eddington, Eddington and sub-Eddington scenarios re- 
spectively. Top Panel: mass function of accreting DCBHs at z = 6 for 
radiative efficiency e = 0.1. Bottom Panel: cumulative mass density of the 
accreting DCBHs vs. redshift for e = 0.07, 0.1, 0.2 represented by dotted, 
solid and dashed lines respectively and following the same colour coding as 
the top plot for / c dd- The black circless represent the total mass density of 
the newly formed DCBHs at each redshift. The solid purple triangle marks 
the claim of Tl 1 at z « 8. 



DCBH-hosting and non-DCBH hosting haloes with their respec- 
tive sourc es. Note that in our case, RR(d) is the same as DR(d), 
as used bv lLi etalj|2012T) . since the position of the halo is arbitrary. 

Following the prescription described above, we first construct 
a cross-correlation function ^Haio(DC) which is computed over all 
the newly formed DCBHs, in a given redshift range, and their re- 
spective sources. We then define a similar cross-correlation func- 
tion CHaio(NoDC) f° r haloes that do not host a DCBH with all their 
respective sourcesQ Both £ H aio(DC) and £ H aio(NoDC), are con- 
structed using the same bins at a given snapshot with the first bin 
placed at a distance larger than the typical virial radius (~ 1 kpc) of 
a massive halo. This is done in order to exclude sources that formed 
within a non DCBH hosting halo at earlier times. This choice does 
not affect the nature or trend of CHaio(NoDC) • 



6 The non DCBH hosting haloes are chosen in the same mass range as the 
haloes hosting a DCBH. This is done for consistency in the construction of 
the cross-correlation function in the two cases. 
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Finally, to check for a variation in the clustering of sources 
around DCBH hosting haloes versus non-DCBH hosting haloes, 
we define 



£Halo(DC) 
£Halo(NoDC) 



(24) 



If CHaio(DC) = £Haio(NoDC) in each distance bin, it implies that the 
DCBH host haloes are as clustered as the non-DCBH host haloes at 
all scales. Therefore, a value of unity for ftotal at a given distance 
scale would imply the lack of clustering for the DCBH host haloes. 
Additionally, a variation in £totai with distance would imply the 
difference in clustering of the DCBH host haloes vs. non-DCBH 
host haloes with the neighbouring sources. A negative slope would 
imply over-clustering for DCBH host haloes at smaller distances 
whereas a positive slope would imply an over-dense environment 
at smaller distances for the non-DCBH hosting haloes. The results 
of our calculations for the cases esc0.5 and esc 1.0 are plotted in 
Fig. [100 

It can be inferred from the value of £totai and the slope of the 
fits that the haloes that host a DCBH are more clustered than the 
non-DCBH hosting haloes especially at a scale of few tens of kpc. 
This trend can be attributed to the fact that in order to reach J cr it 
near a pristine halo, a source (or a population of sources) must exist 
very close by. 

The different escape fractions make this trend even more 
prominent since in both the cases, the line gets flatter as we move 
to lower redshifts but the relative change in the slope of the lines 
is inversely dependent on the escape fraction. This can be under- 
stood by noting that to produce the same level of LW radiation in a 
neighbouring halo, the sources would need to be closer to the halo 
(hence more clustered) if the escape fraction was set to 0.5 instead 
of 1.0. This also implies that the haloes that host DCBH in esc0.5 
form from regions of higher over-densities than in esc 1.0. Hence, 
the flattening trend is most pronounced in the esc 1.0 (red line) as 
compared to the esc0.5 (black line). 

By using our detailed prescription for the spatial variation of 
the LW flux, we find that the haloes which are exposed to J cr i t 
always have a source within a few kpc. We note that our results are 
in accordance with previous work done by D08 which shows that 
the distance scale within which a halo should have a close by LW 
source in order to undergo direct collapse is ~ 10 kpc. 

The function £totai follows a linear fit in log space which can 
be parameterised as 



log(£total) = ci + c 2 log(dphy) 



(25) 



where <f p h y is the physical distance between a halo and LW source 
in parsec and the parameters c\ and C2 are indicated on the bottom 
left of each plot in Fig. 1101 

The flattening of the slope can be attributed to the fact that 
it is more common for a halo to be in an environment with close 
by sources at later redshifts due to the higher overall SFRD in the 
box. Also, the haloes that host DCBHs at z > 10 originate from 
regions of larger overdensities than the ones at lower redshifts. This 
is evident from the larger negative slope of £totai at z > 10. At 
lower redshifts, a lower-a fluctuation is required to produce a halo 
of ~ 10 7 Mq, which roughly corresponds to a T V - 1T = 10 4 K 



7 The chosen cases have at least a few DCBHs in the specified redshift 
bins. The case esc0.5Reion and escO. 1 have very few DCBHs and the case 
0.5HSFE has not been plotted just to avoid repetition as it lies between the 
two plotted cases. 
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Figure 10. The correlation function ftotal plotted for a range of redshifts 
(top right of each plot). Black crosses (line) and red squares (line) indicate 
the correlation data points (fit to the points) computed from the function 
ftotal f° r me cases esc0.5 and escl.O respectively. The numbers on the 
bottom left indicate the fit parameters ci and c 2 in Eq |25| 



halo, supporting that DCBHs must arise from regions of high over 
densities in the early Universe. 

An important inference can be drawn from the points above. 
If the DCBHs which form early on (at z > 10) in the Universe 
arise from highly clustered regions, their environment can only get 
more clustered as the halo progresses to later times. This is an im- 
portant result as it supports the idea that the most massive SMBH 
that we observe at the centres of ellipticals or as highly clustered 
AGNs, might have originated from regions of high over-densities 
quite early on in the Universe, possibly as DCBHs. 



3.4.2 History 

It is interesting to check when the haloes that host DCBHs first 
originated during their cosmic evolution. We checked the most 
massive progenitor of each of the DCBHs and tracked it back in 
time until the halo was found to have a mass equal to the mass res- 
olution in our work (20 DM particles or 1.8 x 10 5 Mq). We label 
this time as the time of the halo's birth, t^ iTth . We define the age of 
the halo when it was first found to host a DCBH as 



,h ±h 

T — l DC ~ I birth 



(26) 
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Figure 11. Age distribution of the DCBH host haloes in our three main 
cases. 



where t^, c is the time when a halo is found to host a newly formed 
DCBH. 

The histograms of r for our 3 main cases are plotted in in 
Fig. [TT] In the case of escO. 1 (blue), 90% of the haloes hosting a 
DCBH were born within 150 Myr, with the remaining haloes being 
500 Myr old. However, for esc0.5 (black) and esc 1.0 (red), all the 
DCBH host haloes are distributed over the age parameter. Part of 
the reason for this changing trend is that a larger / esc implies a 
faster build-up of the LW flux to larger values which leads to a 
higher suppression of Pop III star formation over a longer time. 
The aforementioned suppression of Pop III star formation would 
also be in minihalos formed at earlier times which could later grow 
into massive haloes and form a DCBH if they have a close by Pop 
II source. Thus while in the case of esc0.5, 41% of the haloes are 
less than 150 Myr old, only 11% of similarly aged haloes are seen 
in esc 1.0. 



3.5 Efficiency of DCBH formation 

We plot the redshift evolution of the efficiency of DCBH forma- 
tion, effj^Q in Fig. [12] The efficiency at a redshift is defined as 
the number of newly formed DCBHs (which by definition form in 
pristine massive haloes) divided by number of newly formed mas- 
sive haloes at that redshift. Note that in our model a newly formed 
pristine halo with T V - 1T ^ 10 4 K at a given redshift will immedi- 
ately form either a DCBH (if Jlw Jcrit) or a Pop III star (if 
Jlw < Jcrit). The effj^Q can be expressed as a function of redshift 



(l + *) e 



(27) 



effnc = e i 

where the fit parameters ei and ei are listed in Table|4]for all our 
five cases. Since the efficiency of DCBH formation is a combi- 
nation of DCBH formation rate and the formation rate of newly 
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Figure 12. The efficiency of DCBH , measured as the ratio of the number of 
DCBHs and the total number of newly formed massive halos, as a function 
of redshift. The fit parameters to the lines are listed in Tableg] 



Table 4. Fit parameters to Eq |27l for each case. 



Case 


ei 


e2 


escO. 1 


0.055 


-0.26 


esc0.5 


0.12 


-0.51 


esc 1.0 


170 


-3.32 


esc0.5HSFE 


1.4 


-1.29 


esc0.5Reion 


0.0023 


1.1 



formed massive haloes (which is the same in all the cases), the 
trends in Fig.[T2]are similar to the ones in Fig. [8] Again, the same 
reasoning that applies to the DCBH formation rate, applies to the 
trends in the efficiency. A higher number of LW photons leads to a 
higher efficiency of DCBH formation. The case 1.0 has the highest 
efficiency of DCBH formation followed by esc0.5HSFE which is 
due to the fact that a higher output of LW photons is seen in the 
former case than the latter at z < 11 (see Fig. [5}- However, the 
efficiency in esc0.5Reion decreases at later times which is in accor- 
dance with the low overall DCBH formation rate for this case and 
due to the flattening in the formation rate of pristine haloes we find 
at later times. Note that the similar values of the fit parameters for 
escO. 1 and esc0.5 is due to two reasons: the formation rate densities 
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of DCBHs in these two cases are similar and the formation rate of 
massive pristine haloes in all the cases is exactly the same since it 
is drawn from the same A'-body simulation. 

For the first time, we are able to constrain the seeding mech- 
anism of BHs at high redshifts by self consistently accounting for 
the physical processes that give rise to the conditions for massive 
BH seed formation. Equation [27] encapsulates information about 
the number of massive metal free haloes appearing at a given epoch, 
their clustering with the sources (or haloes) and the rate at which 
these DCBHs form. 



3.6 Reionisation Feedback 

We ran our fiducial case with the ad dition of a simple rei onisation 
feedback prescription motivated by iDiikstra et al.l d2004l) . During 
reionisation, the atomic H ionising photons ionise the gas in mas- 
sive haloes (predominantly comprised of atomic hydrogen) which 
results in the delayed collapse of gas. As a consequence a larger 
potential well is required for gas collapse due to the added gas pres - 
sure from photoheating. In accordance with IDiikstra et all d2004l) , 
a circular velocity threshold of V c = 20 km s _1 was added on top 
of the fiducial model to account for the reionisation feedback in 
haloes with T vir > 10 4 K between 6 < z < 11. In essence, this 
models the impact of instantaneous reionization at z ~ 11. 

The results reflect an immediate quenching of Pop III SF at 
z < 11 in massive haloes in Fig. [4]. The mini-haloes are also unable 
to make Pop III stars due to the high Jb g already in place but the 
Pop II SFR remains almost unaffected. This is because our Pop II 
SF threshold already requires a halo to have A/i n f a ii > 10 8 Mq 
which roughly translates to a V c ~ 20 kms -1 . 

The most interesting outcome is the appearance of only 4 
DCBHs in our box, as compared to the 59 in our fiducial case of 
esc0.5, with no DCBHs seen between 8 < z < 11. This accounts 
for the effect where even though a pristine massive halo in this red- 
shift range might be exposed to J cr it, most of the gas would be 
in a hot ionised state which would prevent it from collapsing and 
forming a DCBH. Only once the halo has a circular velocity greater 
than 20 kms -1 , the gas inside it can collapse and form a DCBH, 
which happens in our box at z < 8. The DCBHs that form before 
the onset of reionisation at 2 > 11 are the ones found in pristine 
massive haloes with no constraints on their circular velocity. This 
is one reason why even though the green and black lines trace each 
other in Fig. [5] only 4 DCBHs are seen in the esc0.5reion case as 
compared to the 59 in the esc0.5 case. 

We allowed the DCBHs in the esc0.5Reion run to grow in the 
same way as the in esc0.5 run (described in Sec. 13. 3t . The early 
appearance of the 3 DCBHs at z > 11 in this case allows them to 
grow into SMBHs (see Fig.l 13b by z ~ 6 with / ct jd = 1. 

We match the recent claim made by Til for a population of 
obscured IMBHs, at z ~ 8, by setting the / c( jd = 0.1 and e = 0.2 
or /edd = 0.4 and e = 0.1. It is interesting to that we are able 
to match Til's claim with a sub-Eddington efficiency in the case 
where we find the least number of DCBHs. 

The fact that DCBH-hosting halos are clustered, also suggests 
that they form in regions of the Universe that are reionised at rel- 
atively early times, as well, due to the concentration of ionising 
sources around them. This suggests that the feedback from reion- 
isation could be even stronger than what we have found assuming 
instantaneous reionisation at z = 11. 
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Figure 13. Same as Fig. [9] but for the esc0.5Reion case. The presence of 
only 4 DCBHs has severe consequences on the BH mass function. We are 
able to match the claim of Til, marked as the solid purple triangle, with 
both Eddington and sub-Eddington accretion modes. The zero points in the 
mass function arise due to the lower (factor of 10) number of DCBHs in the 
esc0.5reion case as compared to esc0.5. 



4 OBSERVABILITY OF THE STELLAR SEEDS OF 
DIRECT COLLAPSE BLACK HOLES 

We have found that a significant number of direct collapse black 
holes are likely to have formed in the early Universe. Now we turn 
to the question of whether these objects are plen tiful enough for 
future surve ys to detect t hem. A s discussed by e.g. lBromm &Loebl 
d2003h and iBegelmanl feoich . the hot protogalactic gas is ex- 
pected to first collapse to a supermassive primordial star which 
subsequently accretes gas until it attains a mass of > 10 4 M 
and collapses to a b l ack hole (see also iDotan & Shavivl 120121 : 
iHosokawa et~aHl2012l ; I Johnson et ai]l2012n . Here we focus on the 
prospects for uncovering these supermassive stellar progenitors of 
direct collapse black holes, as these objects a re expected to be ver y 
bright and possibly detectable by JWST (e.g. lGardner et al.ll200r3) . 
We shall address the question of the detectability of accreting direct 
collapse black holes in future work. 

In order to estimate the likelihood that a given deep survey 
could find SMSs, the precursors of DCBHs, we use the fits provided 
in Table 3 to the rate \dN/dz\ of SMS formation (equal to the rate 
of black hole formation) shown in Fig. 7. With this, we find that the 
expected number nsMS of SMS that lie within a region of the sky, 
as a function of redshift z, is given by 
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dz 
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idJV, 
dz I 



10- 3 Mpc- 3 J Vl0 6 yr 



(28) 



where dV/dz is the comoving volume element per unit redshift, 
\dt/dz\ is the rate of change of the Hubble time with redshift, and 
tiif c is the lifetime of a SMS. To obtain the second equation above 
we have neglected the effect of dark energy on the rate of Hubble 
expansion, which is a reasonable assumption at the high redshifts 
(z > 6) we are considering here; otherwise, we have adopted the 
same cosmological parameters as described in Section 2.1. Finally, 
note that the longer the stellar lifetime tme, the more objects will 
be visible within a given redshift interval. 

Fig- 031 shows the number of SMS per square degree per red- 
shift interval that we find for each of the five cases shown in Fig. 7, 
normalized to tm e = 10 yr, a typical value expected for a rapidly 
accreting SMS (see iBegelmanl |2010| : fjohnson et al]|2012l) . Also 
shown is the minimum number of SMS that would yield an average 
of one SMS per redshift interval (Az = 1) within the area of sky 
cove red by the Deep- Wid e Survey (DWS) planned for the JWST 
(e.g. lGardneretaf1l2006l) . Clearly, the prospects of detection are 
good, as in each of the cases we find that at least a few SMS should 
lie within the survey area ~ 100 arcmin x arcmin. 

We note that rapidly accreting SMS are expected to have 
distinct observational s ignatures which coul d be detected by the 
JWST, as discussed by iJohnson et all J2012h . In particular, these 
objects may emit strong continuum radiation below the Lyman 
limit, and they are likely to exhibit both strong Ha and He II 
A 1640 recombination line emission. An important difference be- 
tween these objects and others with strong recombination line emis- 
sion is that they may also be very weak Lya emitters, due to the 
trapping of Lya photons in the optically thick accretion flows feed- 
ing their growth. The detection of objects exhibiting these observa- 
tional signatures would provide important constraints on both the 
nature and abundance of the stellar seeds of direct collapse black 
holes. 
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Figure 14. The number of supermassive stellar progenitors of direct col- 
lapse black holes, as observed on the sky per square degree per redshift 
interval (Az), as a function of redshift z. The number of supermassive stars 
is given by the fits shown in Fig. 7 and in Table 3, for each of the five cases 
shown. The gray dotted line shows the number of supermassive stars that 
must be present for at least one per redshift interval ( Az = 1) to appear in 
the field of view of the Deep-Wide Survey planned for the JWST. For all 
cases the survey should be large enough for at least a few to of the order of 
10 supermassive stars to lie within the field of view. Strong continuum and 
Ha and He II A1640 emission lines may be detected from these objects. 



exceed J™ t . 

(ii) We find the first DCBH at z « 12, however the total number 
of such objects depends on the LW photon output of a given model. 

(iii) In each of our cases, all the haloes that host DCBHs have 
close by LW sources within ~ 10 kpc. 

(iv) We also find that the haloes that host DCBH at z > 10 are 
more clustered with external LW sources than the haloes that do 
not host DCBH. Also, the DCBHs that appear later at z < 10 are 
less clustered than the DCBHs that appear at z > 10 and may exist 
at the centres of galaxies of various morphological types at z = 0. 



5 SUMMARY AND DISCSUSSION 

In this paper, we present the results from a A^-body, DM only sim- 
ulation of a 3.4 Mpc/i" 1 co-moving box from cosmological ini- 
tial conditions. On top of this simulation, we developed a SAM 
which takes into account the self consistent global build up and lo- 
cal variation of the LW radiation field due to the Pop III and Pop II 
stellar sources. The merging histories of haloes are also tracked in 
order to account for metal pollution form previous episodes of star 
formation. This allowed us to identify the possible sites of DCBH 
and to investigate their environment. Though our simulation is not 
large enough to probe a wide range of environments, we show that 
even in such volumes, BH seeding by DCBHs could be a com- 
mon phenomenon. Our study in this respect motivates the seeding 
of present-day SMBHs via the formation of DCBHs. The key find- 
ings of our work are summarised below. 

(i) DCBH formation sets in with the onset of Pop II star 
formation as the Pop II stars can easily produce Jlw ^ Jllit- 
On the other hand, LW radiation from Pop III stars is not able to 



(v) In our model including reionisation, we are able to match 
recent claims made by Til about the population of obscured 
IMBHs, by assuming both Eddington (/ c( jd = 0.1, e — 0.2) and 
sub-Eddington (/odd = 0.4, e = 0.1) accretion modes for the 
DCBHs. 

(vi) We find that for all our cases, the JWST should be able to 
detect at least a few of the supermassive stellar precursors of these 
DCBHs over a wide range of redshifts (z > 6). 

Our results are subject to limitations due to the modelling ap- 
proach we chose. The halo threshold mass assumed in our work 
for Pop II star formation sets the clock for DCBH formation. In 
our current work, we have set the mass thres hold for Pop II sta r 
formation to 10 s Mq following the work of lMaio et ail d201 ll) . 
however setting it to a lower value would allow for the Pop II stars 
to form earlier in the box. This would lead to an earlier epoch of 
DCBH formation but it is difficult to predict their abundance at 
later times. Note that we also set the mass threshold for structure 
formation (Pop III, Pop II or DCBH) to V c = 20 kms -1 between 
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6 < z < 11, which only further quenches the DCBH formation 
rate. 

Note that our simulated volume is smaller then typical vol- 
umes probed by current observations and does not include sources 
as luminous as the ones detected in the surveys. Thus our predicted 
SFRD should be somewhat lower than the observed ones. Based on 
the observational constraints, as shown in Fig.|3]and|4] we estimate 
that our computed value of Jb g could be lower by a factor of ~ 2 
at a given redshift. A higher level of LW background would make 
it relatively easy to reach J cr i t and would also imply the quenching 
of Pop III star formation in a larger number of minihaloes. Whether 
this would also lead to a higher number of DCBHs is non-trivial to 
predict. 

In principle, SN explosions from neighbouring ste llar popula- 
tions can enrich a pristine halo early on in its lifetime dMaio et alj 
l201lh . This could also reduce the number of DCBHs we find in 
our study, if the metal en richment is high eno ugh to alter the cool- 
ing properties of the gas dOmukai et alfeOoH) . However, it is very 
likely that the metals carried in the SN wind may not be mixed into 
the dense gas at the centre of the halo where DCBH formation oc- 
curs (e.g. (Cen & Riquelmell2008h within the timescale of DCBH 
formation which is ~ 2 — 3 Myr. 

The gas within the haloes identified as DC candidates would 
still need to collapse without fragmentation into a central massive 
object. The study by LN06 explores a mechanism where a Toomre- 
stable gaseous disc in a pristine low spin halo can effectively re- 
distribute its angular momentum, thereby preventing fragmentation 
and eventually forming a DCBH. The aim of our next study is to 
self consistently explore the mechanism suggested by LN06, on top 
of our existing framework, which should in principle greatly reduce 
the number of DCBH host haloes since low spin halo es at such high 
redshifts are quite rare (eg. lDavis & NataraianfeOloT) . 

The accreting discs of both Pop I II remnant BHs and 
DCBHs could also emit LW photon s (e.g. IPelupessv et alj[2007l : 
I Alvarez et all 120091: Ueon et alj|201ll) . where th e emission woul d 
depend on both the BH mass and accretion rat e jGreif etal.ll2008h . 
However, a recent study ( I Johnson et al.ll201ll) has shown that due 
to the low accretion rate of these BHs, their contribution to the LW 
specific intensity can be quite low outside the halo at ~ 1 kpc, i.e. 
the typical virial radius of BH host haloes at high redshifts. Due to 
the uncertainty in the emission characteristics of BH accretion disc, 
we focus on the stellar components to account for the LW radiation 
in our model. Also, the X-ray feedback from accretion discs could 
heat the gas in surrounding halo es, thereby preventing them from 
collap sing and making stars (e.g. lMirabel et al.ll201ll ; lTanaka et alj 
|2012|) . hindering the formation of DCBHs in the neighbouring pris- 
tine massive haloes. We plan to explore the impact of accreting BHs 
on the formation of DCBHs in a future study. 

Since the plausibility of direct collapse in a pristine halo de- 
pends on the number of LW photons reaching it, we find a clear 
degeneracy in the various cases that span the (/ csc ,a) parameter 
space. The degeneracy in the number of LW photons produced in 
our cases could be broken by comparison of the BH mass function 
or the mass density of BHs that we find with those inferred for BHs 
from observations at z > 6. Another possibility is via the detection 
of SMSs in the planned surveys of the JWST which could also shed 
light on the plausibility of this scenario as a SM S is believed to be 
the precursor of a DCBH (e.g. lBegelman]|2010l) . 

Our results shed new light on the long-standing argument that 
only very close star-halo pairs could give rise to DCBHs and that a 
characteristic length of ~ 10 kpc is the maximum distance within 
which a halo must see a LW source in order to have direct col- 



lapse of gas (D08). We find that although a source must exist within 
10 kpc, it is not necessary that a single source produces all of the 
LW radiation which accounts for J cr i t ; there is a contribution from 
the cosmological background LW radiation field, as well as from a 
number of local sources producing J cr it- However, in all the cases, 
it is only the Pop II star clusters that produce all (or most) of the 
J cr it. Pop III stars alone never produce enough LW photons to 
achieve J cr it- 

It is interesting to note that even in the worst-case scenario 
for the formation of DCBHs i.e. the model including reionisation 
feedback (esc0.5Reion), we still find a few DCBHs, which hints to- 
wards the high plausibility of the DCBH scenario. While photoion- 
isation strongly inhibits the formation of Pop III stars in smaller 
pristine haloes, it still allows for Pop II star formation in massive 
enough enriched haloes, which produce the necessary background. 

Allowing the DCBHs to grow via different modes of Edding- 
ton accretion gives rise to a range of possibilities for the BH mass 
function, and can readily account for the presence of supermassive 
black holes by z = 7. The expected number of SMBH is a few 
per co-moving G pc 3 , in accordance w i th the inferred numb er of 
quasars at z > 6 dFan et al ]|200il2006l : lMortlock et al.ll201ll) . We 
over predict the number of such SMBHs but argue that our work is 
an upper limit f or the existence of such objects. However, the study 
by Tl 1 (see also lWillottl201~il : lFiore et al.(2012l) suggests the possi- 
bility of a large number of intermediate mass black holes at z > 6. 
They infer (via extrapolation) the presence of an obscured popu- 
lation of intermediate mass BHs by looking at the stacked X-ray 
luminosity signals of high redshift galaxies. We are able to match 
their claim at z ~ 8 in our reionisation model, assuming both Ed- 
dington and sub-Eddington accretion modes for the DCBHs. Inde- 
pendent of the claim made by Til, on the basis of our model we 
argue that a population of BHs must be present at z > 6 due to the 
sheer number of DCBH host haloes that we find. 

A precise seeding mechanism of BHs at early redshifts in cos- 
mological simulations is important in order to explain the AGN lu- 
minosity functions, growth of massive BHs and the evolution and 
properties of galaxies at lower redshifts. The environment of these 
BHs would play an important role in determining their evolution 
and the Eq.|27]is a first step towards constraining the environments 
and masses of seed BHs. Using a semi-analytical model, which 
takes into account the halo histories and the spatial variation of the 
LW flux, we were able to parameterise the fraction of newly formed 
haloes with T vir > 10 4 K that are able to host DCBH as a func- 
tion of redshift. The equation is an outcome of our model where 
we are able to resolve haloes with masses in the range 10 6 " 7 M 
and could thus serve as a sub-grid model for the seeding of BHs 
in large scale cosmological simulations, which we will pursue in a 
future study. 
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Figure Al. The mass function of haloes at z = 6. The red dashed line de- 
picts the mass function of subhaloes and blue dashed-dotted line represents 
the FoF mass function in our work. 



APPENDIX A: DETAILS OF METHODOLOGY 
Al Mass Function at z — 6 

We plot the mass function of the haloes in our work at z — 6 in 
Fig. IA1I The red and blue lines depict the subfind and FoF halo 
mass function in our work respectively. It is clear from the plot that 
we probe the low mass end of the mass function at z = 6. 



A2 Selection of LW sources 

The stellar populations (both Pop III and Pop II) are the primary 
source of LW radiation at early epochs. However, only certain 
sources can contribute to a LW radiation background at any given 
snapshot. Two important processes govern the fate of these LW 
photons; first, they might get cosmologically redshifted out of the 
LW band while simultaneously, the photons from the bluer end of 
the spectrum enter the LW range and second, these photons can 
get a bsorbed by the neutr al hydrogen present in the early Uni- 
verse. lHaiman et al.l | |2000|) looked at the absorption of these pho- 
tons by the neutral hydrogen present in the Universe. They con- 
cluded that the Ajzlw over which a LW photon can exist is quite 
small since it gets readily absorbed by atomic hydrogen present in 
the un-reionised Universe. The LW band range lies very close and 
even overlaps with transitions that occur in atomic hydrogen, hence 
the mean free path for a LW photon is smaller than the distance it 
can travel before it g ets cosmologically redshifted out of the band 
dHaiman et alj|2000l Fig. 16). Using the relation 



1 + Zr, 



1 + -Zobs 



Z^obs 



(Al) 



one can easily compute the maximum redshift (z max > z bs) at 
which a photon emitted at frequency v t can contribute to the LW 
band, at a given observation redshift (z Q & s ) for a given observation 
frequency (v bs). The upper limit on the lookback redshift (zib) or 
lookback time (in,) can be obtained by setting Vi — 12.1 eV (owing 
to the Lyman-/? line) and v b s = 11.2 eV for any given redshift. 
Hence while calculating the mean LW background for a given red- 
shift, we only count the stars whose photons originate after iib. To 
do this, we define two important parameters for each star/stellar 
population in our study; the time of formation, which refers to the 
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age of the Universe when the star was formed denoted by £f orm and 
the age of the Universe when the star died denoted by 



tc 



tfoTm H~ ^lifc 



(A2) 



The lifetime tuf e of the star depends on the mass of the star and is 
computed using the fits mentioned in Table Q]for Pop III stars (typi- 
cal ftjf e of a 100 Mq star « 2.79 Myr) and lpadovani & Matteuccil 
i 19931) for Pop II star cluster (typical tufe for a Pop II cluster 
weighted by IMF used ~ 10 Gyr). Hence, the selection criteria 
for sources contributing to the Jb g becomes 



^contrib ^ £lb 



(A3) 



In addition to Eq |A3l the selection criteria for stars that can con- 
tribute locally to the LW radiation level also needs to be consid- 
ered. In order to do this, we use a similar approach to KA09,and 
analyse the past light cone of a halo and compare it the world lines 
of the sources. At a given timestep ti, we check if LW photons 
from a source can contribute to the Ji oca i in a halo by comparing 
the physical distance between the source and the halo's position, 
with the time required for the radiation to travel between the birth 
of the source and ti, and the death of the source and ti. In case of 
emission from Pop II stellar clusters, we also calculate their age in 
order to determine when the photons were actually emitted (Fig.[2]l. 

We describe our prescription for selecting stars that are con- 
sidered to contribute spatially to the LW intensity in a halo at time 
ti by writing the conditions 



d s _h < At, if , 

d s -h D It.ic , 



(A4) 
(A5) 



where d s _h is the physical source-halo distance, ti is the age of 
the Universe at snapshot i, D\ t ,a and Di t ,i c represent the physi- 
cal distance light can travel between ti and tform and between ti 
and tcontrib respectively (see Fig. lA2l for more details). Hence if a 
source/star satisfies Eq s . I A3 1 1 A4 1 and |A5l then it contributes locally 
to the LW radiation level and the selection criteria for the sources 
that can contribute to Jb g is given by Eq |A3l 

It is important to note that every halo early on in the Universe 
is expected to be exposed to a minimum level of Jlw given by 
Eq. [12] and [T3] which is an approximation to the mean-background 
level of radiation that is believed to be present everywhere in the 
Universe. Hence, in Eqs.[l2]and[T3]we assume that the SFR density 
in a ~ few Mpc-side box (and hence the comoving density of stars) 
would be the same everywhere in the Universe (see Fig. [A3). 
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Figure A2. Lightcone diagram for the selection criteria of LW sources in 
our work. Red stars indicate the tf orm and maroon filled circles represent 
the t con t r ib for a stellar source. The halo for which the LW intensity is to 
be calculated is placed at i whose past lightcone is marked in dark blue. 
Although all the sources a,b,c,d,e satisfy the Eq. I A3 1 and will contribute to 
Jbg, only the LW photons from stellar sources a and b can make it to the 
halo following Eq . IA4l and lA5 l and will contribute to Ji oca i. 




Figure A3. A halo is (periodically) placed at the centre of one of the sides 
of the simulation box (red) at epoch ('. The halo's past lightcone is denoted 
by the dark blue lines. The light blue line marks the lookback time com- 
puted using I All The actual background Jlw would come from the entire 
Universe (orange box) which can in turn be imagined as a conglomerate of 
smaller simulation boxes (green). The background Jlw is computed using 
the red box but it is assumed that the same mean Jlw would exist through- 
out the Universe (orange box). 
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